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Ambient noise measured in the deep ocean is studied in the context of a search for signals from 
ultra-high energy cosmic ray neutrinos. The spectral shape of the noise at the relevant high fre- 
quencies is found to be very stable for an extensive data set collected over several months from 49 
hydrophones mounted near the bottom of the ocean at ~1600 m depth. The slopes of the ambi- 
ent noise spectra above 15 kHz are found to roll-off faster than the -6 dB/octave seen in Knudsen 
spectra. A model attributing the source to an uniform distribution of surface noise that includes 
frequency-dependent absorption at large depth is found to fit the data well up to 25 kHz. This depth 
dependent model should therefore be used in analysis methods of acoustic neutrino pulse detection 
that require the expected noise spectra. 

PACS numbers: 43.30.Nb, 92.10.Vz, 95.85.Ry 



I. INTRODUCTION 

The study of high frequency acoustic transients in 
water was first proposed by Askaryan [l| as a possible 
method for the detection of the highest energy cosmic 
radiation. In the last few years feasibiUty studies have 
been initiated for acoustic ultra-high-energy neutrino de- 
tectors in large natural bodies of water, ice and salt. This 
technique would present a number of advantages, among 
them the possibility of building very large arrays (thou- 
sands of km^) with sparse hydrophones, owing to the 
large attenuation length (^1 dB/km[2]) of sounds at the 
appropriate frequencies. This is essential because of the 
expected extremely low neutrino flux (<C 1 km~^ yr~^) 
at the energies of interest (E>10^* eV). Because of the 
small interaction cross section of neutrinos, sound has 
to be detected in rather deep media, providing a larger 
probability of interaction per unit area. The acoustic ra- 
diation is produced from the volume expansion caused 
by the heating of the medium where the neutrino inter- 
acts and stops. Theoretical models fsj and experimental 
measurements |4| show that the signature of this kind 
of event is a single bipolar pulse, with energy concen- 
trated around 10 kHz. Substantial data processing of 
pulse shapes and multi-phone correlations are required 
to extract these signals from the ambient noise and from 
transients produced by specific sources. 

It is broadly accepted that underwater ambient noise in 
the range of 1 to 25 kHz is generated at the sea surface Q . 
However, starting with the empirical study of Knudsen et 
al.Q, the extensive data supporting this hypothesis has 
been collected only in shallow water or over a limited fre- 
quency band (or both). The variability of ambient noise 
with depth is also measured at one or a few frequencies 
and does not provide the broad range spectral shape. It is 
unclear as to whether some level of bulk emission (apart 
from molecular agitation) might be required to reproduce 



'Electronic address: ' naokok@stanford.edul 



the noise behavior at high frequencies and large depths. 
As early as 1962, when attempting to model ambient 
noise with wave action. Marsh noted that a depth depen- 
dent feature should be included [^]. Short Q, following 
Urick's treatment of deep water ambient noise, developed 
a mathematical treatment for the depth and frequency 
dependent features. These results have not been veri- 
fied on large data sets and broad frequency bands. The 
effort to study acoustic neutrino detection naturally pro- 
duces very large, wide-band data sets, collected at large 
depths, making an extensive analysis of ambient noise in 
this regime possible. 



II. THE SAUND II EXPERIMENT 

The second phase of the "Study of Acoustic Ultra- 
high energy Neutrino Detection" (SAUND II) employs 
a large hydrophone array in the US Navy's Atlantic 
Undersea Test and Evaluation Center (AUTEC) lo- 
cated at the Tongue of the Ocean (TOTO) in the Ba- 
hamas. SAUND II is currently the largest test for the 
feasibility of acoustic ultra-high-energy neutrino detec- 
tion. This program follows a general study of the ex- 
pected performance Q, and a first experimental phase 
(SAUND I) using seven hydrophones at the same lo- 
cation jlO|. SAUND II uses 49 hydrophones that are 
digitized in the water, with data transmitted to shore 
over optical fibers. The array spans an area of ~ 20 km 
X 50 km with spacing of 3 to 5 km. Hydrophones are 
mounted 5.2 m above the ocean floor, at depths between 
1340 and 1880 m and are omnidirectional with a flat re- 
sponse (<5 dB) from 50 Hz to 15 kHz. Between 15 kHz 
and the system cutoff of 40 kHz there is a directional 
dependence in the response, with up to 22 dB difference 
between the vertical and horizontal directions at 40 kHz. 
The gains of the 49 channels coincide to within 1 dB. 
Analog signals are regenerated on the shore station from 
the digital data (for backward compatibility) and fed to 
the SAUND II data acquisition system that re-digitizes 
them at 156 kHz. Since low frequencies are not relevant 



Sea State 

FIG. 1: Histogram of daily averaged sea state calculated from 
539 noise spectra taken by 49 hydrophones over 11 24-hour 
periods ranging from July to November 2006. The sea state 
bin includes quiet periods calculated to be under the average 
sea state noise level. Error bars indicate the magnitude of 
the statistical fluctuations. 



for SAUND II, a high-pass RC filter is applied to the ana- 
log data with a 3 dB point at 100 Hz. A real time analysis 
program running on seven computers records candidate 
neutrino events as well as other data of interest. This 
includes a power spectral density (PSD) obtained by in- 
tegrating 6.56 ms of data, sampled at 6.4 mus, every 5 s 
for each hydrophone continuously while the SAUND II 
system is on. The results presented here are an analysis 
of this PSD data, taken for the purpose of understanding 
the ambient noise background. By agreement with the 
US Navy, the SAUND II data acquisition system records 
data only when the array is not otherwise used. From 
July 2006 to September 2007, the system has been run- 
ning under these conditions stably for a total integrated 
time of ~150 days. 



III. AMBIENT NOISE DATA 

Ambient noise analyses are performed on a subset of 
the data, consisting of 11 data periods randomly spread 
between July to November 2006. Each data period con- 
sists of 24 hours of continuous data-taking. In order to 
average-out intermittent sources that are of transient na- 
ture, such as ship traffic and fish feeding, all 17280 PSDs 
in a day are averaged to produce one spectrum for each 
hydrophone and every data period. For frequencies up to 
15 kHz, the resulting 539 spectra all follow the expected 
j-5/3 Knudsen shape. Therefore, the 1 kHz - 15 kHz 
region is used to unfold the sea state conditions. The 
ambient noise spectrum at sea state zero, in units of dB 
re /LtPa^/Hz, is approximated from the average spectra 
presented in Urick @ as 

F(/) = 101og(/-5/3)^94_5 



1 2 3 4 5 6 7 8 910 20 3" 40 



S 3D 







b 


















_^ Fit Region 
















average 


\ 






68% 


\ 






1 1 99% 








V"' line 

















1 2 3 4 5 6 7 8 910 20 _ 3a, „ 40 

FreqlkHz] 



FIG. 2: Pss-subtracted power spectral density. The average 
of the daily (5 min) averaged spectra over 11 days (24 hours) 
and 49 (1) hydrophones is shown in a (b), along with the 
envelopes containing 99% and 68% of 539 (287) spectra. 

for the spectral region of 1 kHz to 40 kHz. Since the 
dependence of the overall noise level on sea state does 
not influence the f~^^^ falloff, noise levels for higher 
sea states can be found by adding to Eq. ((!]), Pss — 
301og(ns + 1) where represents the sea state Q. The 
resulting analytical form becomes 

P(/,n,) = 101og(r5/3)+94.5 + P,,(n,) (2) 

with a continuous variable Us in P^s as the only free pa- 
rameter to fit each PSD. Data are then fit to Eq. Q using 
the least square method and equally weighing all frequen- 
cies from 1 kHz to 15 kHz. The discrctized value Us is 
used to produce Fig. [T] The fact that July to Novem- 
ber encompasses hurricane season in the Caribbean, and 
hence volatile sea conditions, explains the relatively high 
daily averaged sea states measured. 

After the fitted sea state quantity Pss is subtracted 
from each spectrum, the resulting 539 full bandwidth 
ambient noise spectra exhibit remarkably similar shapes, 
even beyond the 1 kHz to 15 kHz range where the f~^^^ 
power law no longer describes the slopes. Fig. [2^ shows 
the spread in the data by plotting the envelopes contain- 
ing 99% and 68% of the 539 spectra and their average 
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after "sea state subtraction." It is apparent that while 
sea state conditions change drastically, the underlying 
ambient noise spectral shape, even at high frequencies, 
remains consistent. This is significant considering the 
different locations and seasons included in the analysis. 

In order to check that the ambient noise spectra have 
shapes that are consistent over even shorter time scales, 
Fig. [21d plots the spread of the "sea state subtracted" 
spectra averaged over 5 min intervals measured on Oc- 
tober 2nd 2006 at one hydrophone. Again the spread is 
presented in a manner similar to Fig. [2^, plotting the en- 
velope containing 99% and 68% of the 287 spectra along 
with the average. The smooth features of the spectrum 
are stable. The extremely small spread in these spectra 
disfavors the possibility of intermittent sources contribut- 
ing to daily averaged " sea state subtracted" spectra over 
the broad bands analyzed here. This suggests that if 
noise sources other than those causing surface noise give 
an important contribution, they must be continuous in 
time. Very few sources are identified as having a steady 
all-day impact. Thermal noise is expected to become im- 
portant at higher frequencies than measured here, while 
seismic background of continuous disturbance occurs far 
below 1 kHz Noise produced by marine biology is 

expected to have transient and narrow band characteris- 
tics. One well known exception is snapping shrimp beds, 
known for their consistent sound production underwater 
However, there is no evidence that shrimp beds exist 
at these depths, and the similarity of the spectra mea- 
sured over an area of ~1000 km^ disfavours this possibil- 
ity. Therefore, this analysis proceeds to directly test the 
hypothesis that surface generated noise is the only domi- 
nant broad band source contributing to ambient noise at 
these depths and frequencies. 



IV. CORRELATION WITH WIND 

To verify that the overall noise level due to the sea 
state can indeed be expressed by Eq. [21 measured wind 
speeds are compared to the fitted value P^s measured at 
the hydrophone closest to a shore anemometer station. 
Wind speeds, logged every 10 minutes at a site on Andros 
Island at a height of ~15 m, 21 km away from the ver- 
tical of the hydrophone at a bearing of 110°, are shown 
in Fig. [21 In the Figure the direction of the wind and 
the quantity Pgs averaged every 10 min are also shown. 
The particular data period chosen for this study encom- 
passes a rare occasion with a rather stable wind direction 
(~ 100°) over a long period of time (12-hours). The Pear- 
son product-moment correlation coefficient is computed 
between the first 300 min of Pgs data and an equivalent 
duration of wind data, applying different time offsets be- 
tween the two. The correlation coefficient is shown as 
a function of the time offset in Fig. [H The correlation 
coefficient reaches a maximum of 0.85 at a time offset of 
~80 min. This is consistent with the delay expected from 
the distance, approximate wind speed and direction. In 
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FIG. 3: Wind speed and measured offset from the sea state 
zero, Pss, as functions of time starting July 3rd 2006 at 
0124 UTC (top panel). Wind direction is also plotted in de- 
grees North (bottom panel). Wind measurements are taken 
at a site ~21 km away in the range from the hydrophone used 
here. 



order to understand the statistical significance of this re- 
sult, various 12-hour periods of wind data are randomly 
chosen from the 2006 year, and the same analysis is re- 
peated using the original Pss time series from Fig. [31 For 
each 12-hour period, a plot similar to Fig. [His produced, 
and the maximum correlation coefficient within the to 
420 min offset is chosen. Of 100 such 12-hour wind data 
periods, only once a correlation coefficient higher than 
0.85 was observed. Because the wind direction is usu- 
ally quite variable over the travel time between the hy- 
drophone location and the anemometer, conditions like 
those selected in Fig. [3l are rarely found. In addition 
there is no guarantee that the wind front stays coherent 
over the ~1 hour travel to the anemometer station. De- 
spite these limitations, however, it appears that a clear 
correlation is observed between the quantity Pss and the 
wind strength. 



V. DEPTH DEPENDENT AMBIENT NOISE 
SPECTRUM 

Shown in Fig. [5] is the average of 539 ambient noise 
spectra measured without subtracting the increased noise 
level at different sea states. The expected Knudsen form 
of Eq. ^ is also plotted with a Pss value that fits the 
average best between 1 kHz and 15 kHz. This analyt- 
ical spectrum reproduces the data well in this range. 
Above this frequency, however, it substantially overes- 
timates the noise. Next, a depth-dependent modification 
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FIG. 4: Correlation coefficient between the first 300 minutes 
of noise data from Fig. |3] and the equivalent period of wind 
speed data starting at different time offsets (horizontal axis) . 
The main maximum at a time offset of ~ 80 min is followed 
by substantially less significant maxima due to the spikes of 
the quasi-periodic wind at times between 500 and 600 min. 



developed by Short Q is used to modify Eq. ^ to ac- 
count for the large depth. Short's model considers the 
frequency dependent absorption that becomes significant 
when hydrophones are placed in deep waters, and thus, 
the distance from the surface where the noise originates 
becomes large. This modification is particularly impor- 
tant at high frequencies as absorption becomes substan- 
tial. The effective noise intensity per unit band received 
by an omnidirectional hydrophone located at depth h, is 
given by 



Jo{ah) — 27: J a 



7r/2 



sine'd6' 



(3) 



where J^o is the amplitude of the average intensity per 
unit band per unit solid angle radiated by a unit sur- 
face area, 9 is the angle of the ray arriving at the hy- 
drophone measured from the upward vertical direction, 
and a = a{f) is related to the sound absorption coeffi- 
cient a(/) by ah = — 10 log(e~''''). The index n is 1 (2) 
for surface monopole (dipole) sources. Short calculates 
that the noise field below the horizontal plane of a hy- 
drophone {9 > 7r/2) at a height 91 m above the ocean 
bottom is significantly smaller than that above the hor- 
izontal and can be neglected. Therefore, Eq. [3] is inte- 
grated up to = 7r/2 only taking into account direct 
paths from sources. Other assumptions made in deriving 
Eq.[3]are that the noise at the hydrophone is the incoher- 
ent sum of all intensities arriving, after attenuation, from 
the surface sources and that straight ray propagation is 
adequate ^Sj]. Thus, the depth dependent correction to 
the Knudsen spectra can be determined by adding the 
following expression to Eq. ([2]). 



In evaluating this frequency dependent offset, the average 
depth of the 49 hydrophones, h — 1631 m, is used. The 
sound absorption in sea water is parameterized according 
to Fisher and Simmons to evaluate a{f) and therefore 
a{f). A temperature profile taken at TOTO every 7.6 m 
down to 1830 m is used in this evaluation. The resulting 
offset is always negative, as expected for an attenuation, 
and is applied to Eq. The result is also plotted in 
Fig. El 
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FIG. 5: The average ambient noise measured, plotted in log 
(above) and linear (below) frequency scale with 3 theoretical 
curves of Knudsen spectra, as discussed in text. The theo- 
retical curve for thermal noise 0] is also plotted. The curves 
labeled "Short corrections" in the legend include frequency 
dependent attenuation parametrized by Fisher and Simmons 



An attempt to further improve the theoretical curve is 
made by considering the frequency-dependent direction- 
ality of the hydrophones above 15 kHz. This is obtained 
by modifying Eq. [3] into 



.7r/2 

Jeff (ah) = 27rJoo / ' 
Jo 



10 log[ Jo (a/i)/ Jo (0)] 



(4) 



(5) 

with the response function g{6, f) provided as a look-up 
table by the hydrophone specifications. Eq. 2] is then 
evaluated numerically and plotted in Fig. [51 It should 
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be noted that the Pss fit to the data in both of the 
modified Knudsen spectra is performed after the Eq. 3] 
correction and therefore, the numerical values of Pss are 
slightly different in all 3 curves. Also noteworthy is the 
fact that although using n — 1 and n = 2 create very 
small differences as also shown by Short the first case 
(monopole) produces a slightly better fit to data and is, 
therefore, used in the figure. However, the quality of data 
is not sufficient to significantly distinguish between the 
two models. 

Clearly the curves including attenuation provide better 
descriptions of the data at high frequency. Above 25 kHz 
the data appear to flatten faster than the model. It is 
unclear as to whether this should be considered a real 
underwater effect or simply due to difficulties in modeling 
the response of the hydrophones. It should be pointed 
out that the levels are too high to be considered as onset 
of thermal molecular noise, also shown in Fig. [5l 

VI. CONCLUSION 

The analysis of a large acoustic noise dataset collected 
at depths in excess of 1500 m and frequencies up to 
40 kHz confirms that, even at substantial depths, a uni- 
form surface distribution of sources alone is sufficient to 



describe ambient noise characteristics well up to 25 kHz. 
A correlation observed between wind speed and overall 
noise levels at these depths during steady wind condi- 
tions further confirms this finding. A description of the 
spectral shape above 15 kHz is found to require, at the 
depths considered here, proper account of the sound at- 
tenuation in water. This depth effect must therefore be 
included when defining parameters in acoustic neutrino 
detection, such as the transfer function and the signal to 
noise ratio, that require the expected noise spectrum. 
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